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I derive the distribution of energy and momentum transmitted from a fast parton to a medium of 
thermalized quarks, or the source term, in perturbative thermal field theory directly from the quark 
energy-momentum tensor. The fast parton is coupled to the medium by adding an interaction 
term to the Lagrangian. The thermal expectation value of the energy-momentum tensor source 
term is then evaluated using standard Feynman rules at finite temperature. It is found that local 
excitations, which are important for exciting an observable Mach cone structure, fall sharply as a 
function of the energy of the fast parton. This may have implications for the trigger pr dependence 
of measurements of azimuthal dihadron particle correlations in heavy-ion collisions. In particular, 
a conical emission pattern would be less likely to be observed for increasing trigger pr- I show 
that the results presented in this paper can be generalized to more realistic modeling of fast parton 
propagation, such as through a time dependent interaction term, in future studies. 

PACS numbers: 12.38.Mh,25.75Ld,25.75.Bh 



(N 



■ 

Oh 1 
i ' 

43 



(N 
> 
ON 

On 



X 



At sufficiently high temperature, quantum- 
chromo dynamic (QCD) matter is expected to undergo a 
transition from colorless hadrons to a state of deconfincd 
quarks and gluons known as the quark-gluon plasma 
(QGP) jjj. Experimental results from the Relativistic 
Heavy-Ion Collider (RHIC) indicate that the QGP 
has been formed in relativistic heavy- ion collisions Q. 
The results suggest that the QGP formed at RHIC 
center-of-mass energies may behave as a nearly ideal 
fluid Another striking result [|[ is that highly 

energetic, or fast, partons appear to lose a significant 
amount of energy to the medium in a process known 
as jet quenching [5(. In the light of these observations, 
an interesting problem is to calculate how the QGP 
responds to a propagating fast parton. 

This problem has gained attention due to experimental 
measurements of azimuthal particle correlations associ- 
ated with high pt triggers in heavy-ion collisions that dis- 
play a double-peaked or conical structure @ . These mea- 
surements may reflect the interaction of the jet with the 
medium and proposed explanations for the structure have 
included large angle gluon radiation Q, Cerenkov radia- 
tion, Q , colored wakes @ , and perhaps most commonly, 
Mach cone Shockwaves excited in the bulk medium by 
fast partons [l(],[ll|. Other explanations, which do not 
reflect the interaction of the jet with the medium, such 
as fluctuating initial conditions and triangular flow have 
also been proposed [n| . Whatever the nature of the con- 
ical structure associated with high px triggers at RHIC 
is, the QGP response to a fast parton remains an im- 
portant topic, particularly in the light of the surplus of 
high transverse momentum probes available to the newly 
online heavy-ion program at the Large Hadron Collider 
as compared to RHIC. For this reason, the experimental 
handle on medium response to fast partons should only 
improve. 

Theoretical investigation of the QGP response to fast 
partons has mostly been phenomcnological, with many 



studies using hydrodynamics to model the medium re- 
sponse [l3l |. Even within the framework of hydrody- 
namics one still must specify how the fast parton excites 
the medium, that is, the source term. The source term, 
here denoted J v , couples to the energy-momentum ten- 
sor (EMT, denoted T» v in what follows) as d^T^ = J v 
and describes the flow of energy and momentum density 
between an external source and the medium. Previous 
studies have found that the azimuthal dihadron corre- 
lation spectrum associated with a fast parton depends 
sensitively on the form of the source term [To|, [HI • ^ has 
been shown that diffusive momentum, which flows in the 
direction of the source parton's velocity, tends to fill up 
any double-peaked conical structure in the final spectrum 
15, 16]. These observations suggest that the shape of az- 
imuthal dihadron correlations may be a powerful probe 
of the microscopic interaction of energetic partons with 
the QGP. There has been little rigorous theoretical in- 
vestigation of the source term generated by a fast parton 
in QCD. However, a calculation of the source term gen- 
erated by a fast parton for a perturbative QGP within 
kinetic theory was performed in [I?} ■ 

Ideally one would like to calculate the QGP response 
to a fast parton from first principles, such as in the inves- 
tigation of heavy quark propagation through a strongly 
coupled thermal plasma done by Friess et. al using the 
AdS/CFT correspondence [HI- The authors evaluated 
the stress-tensor within the context of linearized gravity. 
It was later shown that the solution matched well with 
a nearly ideal linearized hydrodynamic calculation up to 
distances of about 2/T away from the heavy quark [l9(. 

Performing a similar study within QCD currently 
seems out reach. The problem is inherently dynamic and 
must be approached with real-time (and not imaginary- 
time) methods. It is thus is not well suited for lattice 
techniques, which are more effective in imaginary time. 
Even for a truly perturbative QGP, it is not clear how 
to sum all of the diagrams necessary to calculate the 
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FIG. 1: Feynman diagrams contributing to {d^T^ (x)) in the presence of a source interaction term, A a ^j^. The diagram in 
Figure 1(a) can be traced back to terms in the energy- momentum tensor (see equation fT}) which go as ip-fdip, whereas the 
diagram in Figure 1(b) originates from terms that go as g^Atp. The contribution from the energy-momentum tensor and 
external current insertions are discussed in the text (see also Figure [2j. 



medium response in regions far from the fast parton. In 
view of these limitations, one can consider two (of per- 
haps several) possibilities: a) calculate components of the 
EMT within perturbation theory to obtain information 
about the medium response in regions near the fast par- 
ton; b) calculate the source term for the EMT within 
perturbation theory and use an effective theory to prop- 
agate the resulting disturbance to regions far from the 
fast parton. For this second option, one would calculate 
d^T^ = J" using the fundamental EMT of QCD, and 
then allow d^T^ = J v , where is the hydrodynamic 
EMT, for instance, and then solve the resulting equations 
of motion. 

Of the two options presented above, a) is under more 
theoretical control and will be considered in a future 
work. In this work, however, I consider option b) be- 
cause it is a more straightforward calculation and makes 
more direct contact with previous studies. Furthermore, 
I will consider a medium of massless quarks/antiquarks, 
leaving the inclusion of medium gluons for a forthcom- 
ing study (thus, apart from the coupling strengths, the 
results will be the same as for a QED plasma, for which 
they are also new). The EMT is given by [2(| 

= 1$ (V D v + Y V> - 9^-^, (!) 

where 

JSf=^5v,^ /1 = ^-i5^t° ( 2 ) 

and 

D v ip = D" i> -^7 M D* v ip. ( 3 ) 

In the above equations, g is the strong coupling, t a are the 
SU(3) generators in the fundamental representation, and 
conventional slashed notation is used, 4 = In A 11 , etc. A 
summation over color, spin, and the active number of 
quark flavors is implied in the EMT. 

In order to investigate the medium response it is nec- 
essary to specify how the fast parton couples to the 



medium. One possibility, which will be adopted here, 
is to model the fast parton as an external current which 
couples to the Lagrangian: 

^-►JSf-AJtf, (4) 

where for the moment I do not specify the explicit form 
OI i should not be confused with the source term, 
J u ). The replacement made in (j4]) preserves non-Abelian 
gauge symmetry as long as D^ b j£ = 0. 

The lowest order Feynman diagrams for calculating 
the thermal expectation (d^T^ u (x)) in the presence of 
the interaction term A^jg, are shown in Figure [TJ The 
two gluon exchange is necessary to couple to the EMT, 
which is a color singlet quantity (a two photon exchange 
is also necessary in QED, from Furry's theorem). The 
diagram in Figure 1(a) arises from terms in ([T]) which go 
as ip^dip, whereas Figure 1(b) arises from terms that go 
as g t/jjAifj. My convention is that the standard Feynman 
quark-gluon vertex contributes ig^f^t a , and I will use the 
Feynman gauge for gluon propagators. 

In order to assign a value to these diagrams, one must 
determine what the correct Feynman rules for the EMT 
and external current are. These nonstandard contribu- 
tions are isolated in Figure O The contribution from the 
EMT shown in Figure [2ja) can be determined by simply 
assigning the appropriate momentum to each derivative. 
Recalling that one is here interested in (d^T^) and that 
the final result is Fourier transformed into position space, 
I find that the value of Figure [2j a) is 

£ e — tx-(pi-p 2 ) 

~ (5) 
x ((Pi ~ P i)7" + Pi (3^ 2 +f l )-V2^t l +t 2 ))- 

One can apply the same procedure to Figure [2jb) where 
I will choose the convention that the gluon momentum 
flows away from the external current (or into x and any 
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FIG. 2: Contributions to the Feynman diagrams of Figure 
Q] which come from inserting the energy-momentum tensor 
and external current interaction term. Their values and how 
they are obtained are discussed in the text, see specifically 
equations © - ©. 



vertex) in all cases. The result is 



-ige- i ^ + ^-^(p 4 +p 3 - P2 )^ 



(6) 



which is only valid when the gluon field in Figure Gib) 
connects with the source. 

Finally, for the source contribution in Figure [2lc), one 
has very generally 



(7) 



for the case of an external current which contains one 
power of g. For the sake of simplicity and to be able 
to more easily compare with previous results, I will here 
consider an asymptotically propagating fast parton rep- 
resented by j£ — gQ a (t)U^ <5 3 (z — ut) where u is the fast 
parton's velocity and = (1, u). gQ a (t) is the charge of 
a classical particle in QCD (defined by Q| = 5ijC 2 i, 
with C*2i the quadratic Casimir in representation i (3 for 
a gluon, 4/3 for a quark)) and Q a (t) evolves in time ac- 
cording to Wong's equations [2l|. I will not consider the 



time dependence of Q a {t) here, as it is of higher order 
in g (however, for higher order calculations one must ac- 
count for the time dependence of Q a (t) in order to pre- 
serve gauge invariance). For this choice of j£ Figure Gib) 
simplifies to 



-2irigQ a U a 5{ Pi ■ U). 



(8) 



Before moving on, I point out once again that the Feyn- 
man gauge has been used in obtaining equations ([BE]). 

It is worthwhile to mention here that the above results 
can be modified in a straightforward manner to consider 
more general situations. For example, one could consider 
a fast parton created at an initial time, t = 0, by enforc- 
ing a step function in the external current or even include 
medium induced gluon radiation from the fast parton in 
a phenomcnological way by treating the time dependent 
spectrum of energetic gluons as additional external cur- 
rents in the medium. In any of these situations, the rule 
expressed in (J7J) can be applied. 

As mentioned previously, it is necessary to here work 
in the real-time formalism of finite temperature field the- 
ory. The Green's function structure of Figure [T] is ob- 
tained using the conventions of real-time thermal field 
theory as outlined by Das [22J. In the real-time formal- 
ism at finite temperature one must sum over two types 
of vertices which differ in value by an overall sign. Ad- 
ditionally, propagators separate into vacuum and finite 
temperature contributions. I will not present the details 
of the calculation in this brief report, reserving them in- 
stead for a follow-up publication in which medium gluons 
will also be considered. Using standard Feynman rules 
for real-time thermal field theory in addition to the spe- 
cial rules presented in equations © - ([7]) the result for 
the thermal contribution is obtained as 



(d fi T» v (x)) = -AiN F C 2 g 4 



^^^V ^*^^ . u)s{p5 . v) 



(2tt) 9 



[2 (p 3 • Ufpl - U 2 P3 ■ ViVl - U v ( P3 ■ U)(2 P3 ■ P5 + PA ■ p B )] 
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(9) 



where Gn(p) = (p 2 + iep°) 1 is the retarded Green's 
function, np(p) = (e' p '/ T + l) _1 is the Fermi distribution 
function, T is the temperature, and Np is the number of 
active flavors. 

Equation ([9]) is the central result of this paper and it is 
worthwhile to consider the Green's function structure it 
contains, because this will provide an intuitive picture of 
the underlying physics. If one reads from left to right in 
the first line it is clear that a particle from the heat bath 
(represented by the thermal distribution 5{p 3 )np(p 3 )) 
absorbs a gluon from the external current (represented 



by Gn(pi)) and then continues to propagate (shown by 
Gr{p 3 +Pa)) until it absorbs another gluon from the ex- 
ternal current (represented by Gij(ps)). A final propaga- 
tor representing the reabsorption of the particle by the 
heat bath (which would be given by Gn{pz+Pi+Pb)) has 
been eliminated by the momentum structure of (d^T^). 

In order to further extract meaningful information 
from equation ^ it is useful to consider certain approx- 
imations and limits. As was mentioned above, a kinetic 
theory calculation of the source term generated by an 
asymptotically propagating fast parton was performed 
in [17| ■ Equation © should reduce to the kinetic theory 
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result within the hard thermal loop (HTL) approxima- 
tion (the HTL approximation is formally equivalent to 
the Vlasov equation [23| ). It is worthwhile to check that 
this is indeed the case, as it gives confidence to the ap- 
proach introduced in this paper and that (|9]) has been 
evaluated correctly. 

The HTL result can be obtained in the limit that the 
fields generated by the external current are soft compared 



to the temperature, or formally \p^\ ~ gT <C T in equa- 
tion © (<? is formally considered to be much less than 1). 
Explicitly, one should expand Gr(p3 + Pa) in the limit of 
\p±\ <C |^»3 1 (since \p%\i which appears in the Fermi dis- 
tribution function, is cut off in the integration at a value 
on the order of T) to lowest contributing order. With 
this guidance, I find that in the HTL approximation © 
reduces to 



i m 2 D C 2 g 2 



d pi d f d£l 



2 J (2tt) 6 J 4tt 
pl{{vU) 2 pl~U 1 ' (vU)(v P5 )) 
(p° - v • P4 + it) 2 



-™-^+P*)G R {p A )G R (p 5 )5(p A ■ U)6(p 5 ■ U) 
U 2 vp iP l + U» (vU)( Pi - P5 ) 



VI - V • P4 



(10) 



where m 2 D = g 2 T 2 Np/6 (quarks only). In equation (TTU|) 
dQ is the integration measure over the solid angle de- 
fined by the unit vector v, and = (l,v). Addition- 
ally, bare propagators are used for the Green's functions 
Gr{pa)Gr{p5) in (|10l) in order to facilitate comparison 
to previous results. 

The interested reader can verify (as this author has 
done) that (fTQ|) indeed reproduces the result of fl7| |. 
Rather than prove this for the full spatial distribution, I 
will here consider the total rate of energy transferred to 
the medium and show that this quantity matches with 
the previous result (leading log HTL). The total rate of 
energy transfer, here denoted dE/dt, is found by inte- 
grating the zero component of the source term over all 
space 

— = J d 3 x(d^°(x)). (11) 

The spatial integration yields a factor (27r) 3 <5(p4 + pg) 
from which one can trivially perform the ps integration. 
It is then straightforward to check that the two terms 
with coefficient U v cancel each other, and furthermore 
the term with coefficient U 2 will vanish by symmetry 
upon integration. One finds the expression depends loga- 
rithmically on |p4 1 . The final result including the angular 
integrations is 

dE _ m 2 D C 2 a s A _ tanh -1 [M] \ ^ p max ^ 

dt HTL 2 \ 7 2 U ) p min 

where 7 = (1 — u 2 ) -1 / 2 and based on the approximation 
used in obtaining (|10[) and physical reasoning p max ~ 
T and p m in ~ mo. The expression (|12|) matches the 
leading log result of [l7|. In what follows I will only 
consider the relativistic limit (7 3> 1) of (fT2"|) . 

Having verified that (jlpp reproduces known results 
from the HTL approximation, one can consider the full 
result, ©, which includes contributions from the full mo- 
mentum spectrum (that is, \p^\ ~ \p^\ and \p^\ ^> \p3\). 



In order to simplify the discussion I will consider the 
source term ansatz discussed in [l4| 

(dpT^ix)) ~^§(U"- \d v ) 5(x - ut). (13) 

The utility of (|13|) is that one can encode basic features of 
the source term in a compact way. It is straightforward 
to check that (fTT|) is satisfied by (fl~3|) . The dimension- 
full coefficient A parameterizes local contributions from 
the source, that is, terms which globally integrate to zero 
(and for instance do not contribute to dE/dt) but may 
still be important for exciting the medium. It was found 
in [l4| that within linearized hydrodynamics a double 
peaked structure in the azimuthal emission spectrum as- 
sociated with the source given by (|13|) only appeared for 
rather large values of A (on the order of 0.5 fm or higher 
for 20 GeV total energy deposited into the medium) . The 
basic assumption of equation (|13|) is that the source term 
(|9]) can be expanded in terms of a 8 function and deriva- 
tives of a delta function, which should be a reasonable 
approximation within the context of hydrodynamics. 

It is straightforward to verify that the full source term 
given by equation (|9]) can be put into the form of (fT3|) by 
using (fTTj) to obtain dE/dt and 

A«^ fd*xx(d^{x)) (14) 

where I have chosen to use the x component of the source 
to obtain A (according to the ansatz of (IT51) any compo- 
nent would work) . If one inserts equation (|14l) into (fl3|) it 
is clear that dE / dt drops out of the term proportional to 
the gradient, thus dE/dt is simply a multiplicative factor 
in the evaluation of A (it is written in this way to make 
easier comparison to [1J])- For this reason, and to sim- 
plify the analysis, I will take dE/dt w dE/dtwvL (I here 
remind the reader I use the relativistic limit of dE/dt so 
that it contains no dependence upon the velocity of the 
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FIG. 3: In the source term of equation (|13[) A couples directly 
to sound modes and is important for generating observable 
Mach cone signals. The contribution to A from the source 
term © is plotted above for g = 1 and T = 400 MeV. It 
is clear that A drops as a function of 7 = E/M (M being 
the mass of the fast parton) which could have implications 
for the trigger pr dependence of azimuthal dihadron correla- 
tion measurements (see discussion in text). The inset shows 
that the result for A is largely independent of the separation 
parameter, q* (see text for details). 



fast parton) and instead focus on evaluating the integral 
term in (|14p , for which I will use the full expression given 
by (HI). 

The evaluation of the integral in ([T4"|) can be done di- 
rectly from ^ through an integration by parts of the 
general form: 



J d 3 xd 3 pxe tp x f{p) = i J d 3 xd 3 pe lpx d p J{p). 



(15) 



One has to be careful in using (fT5j) because A has an 
infrared divergence. The simplest way around this, and 
what I will adopt here, is to introduce a mass term in 
the Green's functions Gr(p4)Gr(p5) in © given by the 
thermal gluon mass m 2 = m 2 D /2. 

One can only apply the full result of equation ^ down 
to some momentum scale |j>4 ~ q* where gT « 9' C T. 
When \pi\ ~ gT the full result contains contributions 
from higher order terms and only the HTL approximated 
result of equation (|TU|) is consistent to the order I am 
working. From a practical point of view in the calculation 
of A this means dividing terms in the integration such 
that 



A = Af when \p±\ > q*, 
A = A H tl when |p 4 | < q* , 



(16) 



where Ap is obtained from Q and Ahtl is obtained from 
(fit)]) . The final result should be independent of q* (this 
type of analysis for the separation of scales in calculations 
of finite temperature field theory was introduced in f24j). 



As mentioned above it is not possible to present the 
details of the calculation in this brief report and I will 
instead skip to the result for A and show that it is inde- 
pendent of q* . The result is shown in Figure [3] as a func- 
tion of 7 = (1 — it 2 ) -1 / 2 , and, at the risk of pushing the 
formal limit of g <C 1, I have chosen the experimentally 
relevant values of g — 1 and T = 400 MeV. It is inter- 
esting to note that A falls rather sharply as a function of 
7. Recall that A parameterizes local contributions from 
the source and is important for exciting the medium. In 
particular, A couples directly to sound modes (and not 
to diffusive modes) when (flU|) is used as a source term 
for hydrodynamics and was found to be crucial to the 
appearance of a conical Mach-likc emission spectrum in 
[14j |. The experimental implication of the dependence of 
A on 7 as shown in Figure [3J could be found in the trigger 
Pt dependence of measurements of azimuthal dihadron 
particle correlations. In particular, a conical emission 
pattern would be less likely to be observed for increasing 
trigger pt, which indeed seems to be the case [Hj]. The 
inset of Figure [3] shows that the result is largely indepen- 
dent of the separation parameter, q* , introduced above. 

In summary I have calculated the source term for a 
medium of thermalized quarks in the presence of an 
asymptotically propagating fast parton directly from the 
EMT to leading order in perturbation theory and have 
presented the result in equation ((HJ) . Fcynman rules spe- 
cific to this calculation were presented in equations (|5|) - 
([7]) which can be applied in a straightforward way to more 
complicated external currents in future studies. Within 
the HTL approximation, it was shown that (|9]) reduces to 
a previous result 17| derived from kinetic theory. In ad- 
dition to the HTL approximation, I considered the full re- 
sult, which includes contributions when the source fields 
are hard (that is, \p\ ~ T and \p\ T). In order to 
simplify the analysis, the source was put into the form 
of equation (fT3|) . in which case the source is described 
by the total rate of energy transferred to the medium, 
dE/dt, and local excitations described by the coefficient 
A. I here focused on the evaluation of A and found that 
it falls sharply as a function of 7, as shown in Figure 
[3J This may have implications for the trigger px depen- 
dence of measurements of azimuthal dihadron particle 
correlations in heavy-ion collisions. In particular, a con- 
ical emission pattern would be less likely to be observed 
for increasing trigger pt, which may indeed be the case 
[2o| . A future publication will include medium gluons 
and present the details of the calculation (2(| . 
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